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We study the equilibration behaviour of a quantum particle in a one-dimensional box, with respect 
to a coarse grained position measurement (whether it lies in a certain spatial window or not). We 
show that equilibration in this context indeed takes place and does so very rapidly, in a time 
comparable to the time for the initial wave packet to reach the edges of the box. We also show 
that, for this situation, the equilibration behaviour is relatively insensitive to the precise choice of 
position measurements or initial condition. 


Quantum systems are time-reversible and, in their time 
evolution, have recurrences arbitrarily close to their ini¬ 
tial state. This is in sharp contrast with the irreversibility 
found in thermodynamic processes and with the concept 
of evolution towards equilibrium, which is impossible in 
a periodic system. Lately however, it has been found 
that even a closed quantum system can appear to equili¬ 
brate HHID]. 

Under the weak assumption of a Hamiltonian with 
non-degenerate energy gaps, for a system with high- 
dimension, Linden et al mm proved that the partial 
state of any small subsystem (no matter how you parti¬ 
tion the Hilbert space) stays very close to a static state 
for the vast majority of its time evolution. 

In parallel, Reimann nm showed that, for the same 
set of Hamiltonians, the expectation values of any ‘rea¬ 
sonable’ quantum observable on a closed system also 
stays predominantly close to a static value, even though 
the state itself does not. 

Building on this. Short showed that these results ap¬ 
ply even if one considers the specific possible outcomes of 
the measurement, instead of just the expectation value, 
by describing it as a positive operator-valued measure 
(POVM). Thus, the time-evolving state of a closed sys¬ 
tem is physically nearly indistinguishable from a static 
state, even under careful experimental scrutiny, for the 
majority of time. 

Since then, more effort has been devoted to studying 
the dynamics of the equilibration process |6HI9]. Re¬ 
sults by Short and Farrelly jB] place an upper bound 
on a state’s distinguishability which tightens with time. 
Masanes et al [7] derive the equilibration time scale for 
spin systems averaged over a set of Hamiltonians. Hut- 
ter and Wehner [S] find a clear criterion for when systems 
lose their initial information due to interaction with an 
environment. And Malabarba et al [9] and Goldstein et 
al HD] show that closed systems equilibrate surprisingly 
fast for typical measurements (measurements based on 
random projectors on the Hilbert space). 

Here, we aim to gain further physical insight on the 
precise dynamics involved in equilibration. Showing how 
physical time scales—as we observe in thermodynamic 
systems—emerge from the underlying quantum dynam¬ 
ics has proven a challenge. In fact, it has been shown that 


even for general nanoscale systems you can find measure¬ 
ments where this time scale ranges from nanoseconds to 
the age of the universe laiin]. Therefore, it is clear that 
to solve this problem one must consider scenarios where 
both the Hamiltonian and the measurement are physi¬ 
cally meaningful. 

For this reason, we consider the equilibration of a sim¬ 
ple quantum system composed of a particle in a one¬ 
dimensional box, as perceived by coarse grained position 
measurement (whether it lies in a certain spatial window 
or not). Under these conditions, we study the equilibra¬ 
tion profile and derive the time scales by solving the time 
evolution. 

A special characteristic of the particle in a box is that 
the recurrence time (for any initial condition) is at most 
the ground state period. This allows us to explore the 
entire time evolution numerically, which would be unfea¬ 
sible for systems with very large recurrence times. These 
results are then compared to those cited above and it is 
shown that the equilibration time scale of physical mea¬ 
surements is fast, albeit slower than generic typical mea¬ 
surements. 


I. DEFINITIONS 

The specific system studied is the quantum particle 
inside a one-dimensional box of size L > 0, defined by 
the free-particle Schrodinger equation 


where p is the momentum operator, and the boundary 
conditions restrict the state to the domain x G [^, 
in position space, 

=(x|vk)|^^. =0. (2) 

We consider here a class of projective measurements 
Mw, with w G [0,L]. Each measurement ilAu = 
{Ay^,Byj\ is a set of two projectors in position space, 
each corresponding to a different outcome, defined by 

A^= |a:)(x| dx, Hu, = 1 - Ay,. (3) 
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They represent whether or not the particle is inside a 
region of width w centered on the origin. 

Given a system of dimension d, it can appear to eqnili- 
brate with respect to a particular measurement because 
all the information contained in the system’s degrees 
of freedom (up to d complex numbers) is reduced to 
the probabilities of each measurement outcome (one real 
number for each possible outcome). Under this amount of 
ignorance, even orthogonal states may be indistiguishable 
and the system may appear to be stationary. Choosing 
a two-outcome measurement ensures we are taking full 
leverage of this and taking an approach that is as simple 
as possible to understanding equilibration. 

Our prime object of study is the distinguishability be¬ 
tween two states according to a measurement M 

-Dm(pi,P 2 ) = ^ ^ |Tr[A(pi-p 2 )]| , (4) 

AeM 

which, for reduces to 


II. RESULTS 

First, we consider the initial condition with an approx¬ 
imately Gaussian distribution, 

(x|'kG(0)) = fV ^ 


pG = |^ g )( 4 ' g |, 

where the arbitrary constant cr defines the initial width 
of the wave packet in position space, N is the normal¬ 
ization constant, and the approximate relation is due to 
neglecting the second exponential on the first line. This 
is a smooth function concentrated mostly inside the sup¬ 
port of a single outcome of M {A^), as long as a < w. 
For this system, the eigenstates, the energy levels, and 
the energy amplitudes are 


Dm^{Pi,P 2) = \Ar[Ayj{pi-p2)]\. (5) 

This quantity is significant because it determines the 
probability of successfully determining whether the sys¬ 
tem was originally in state pi or p 2 after performing 
the measurement M |^. This probability is p = ^ + 
\Dm{pi, P 2 ), assuming one has full knowledge of the two 
states being compared. 

The verification that equilibration has happened (at 
some time t) is done by comparing the instantaneous 
state of a system p{t) with its time-average state u) = 
{p{t))te{o 00)1 called the equilibrium state, where 


{x\n) 


(TgIu) 



h‘^'n?Tr‘^ 


2mL2 




( 10 ) 

( 11 ) 


Before describing the distinguishability we state our 
first result about this system. The effective dimension, 
defined in eq. ([^, can be written as 


(/(O)tg(o.oo) = ^ ^ /(Odi. (6) 

Equilibration is achieved once Dm{oj, p{t)) becomes small 
(<C 1) and stays small for the majority of the time evo¬ 
lution. 

Another quantity which plays an important role in the 
equilibration process is the effective dimension EiEiin]. 
Roughly, it represents the number of energy eigenstates 
(|n)) which the system occupies with significant proba¬ 
bility and, for a pure state l^*), it can be written as 


deff 


C30 


-1 




(7) 


Throughout this paper, whenever we mention the high¬ 
dimensional limit we are referring to deff 3> 1. 

Our results, presented in the next Section, consist of 
mathematical statements regarding important character¬ 
istics of the equilibration process of the particle in a box 
system with respect to the simple dichotomic measure¬ 
ment described in eq. ([^. We are effectively asking the 
question: “If all I know is which interval the particle is 
in, does it look like it equilibrates? How fast does that 
happen 



Where in the third line we have used the fact that only 
the odd-n terms contribute to the sum. Furthermore, 
since the error involved in the approximation of eq. ([^ 
is exponentially small in L/ct, for ct/L small enough we 

may incorporate this into O which leads to 

This means that des 1 corresponds to cr ^ L. 
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FIG. 1. Probability distribution of the equilibrium state, 
(xIcugI*) from eq. (141, of the Gaussian initial condition eq. Q 
for several values of deff- 


FIG. 2. Distinguishability, as defined in eq. ||^, between the 
time evolution of pc (eq. ®) and its equilibrium state u)g as 
a function of time, for three different deft and w = ^. Tg is 


the ground state period, see eq. (17l. 


A. The Distinguishability 

The equilibrium state II10E] is calculated as the av¬ 
erage over all time of the initial condition. This is equal 
to the initial state decohered in the energy basis, which, 
for p(0) = Pg reduces to 



where, again, the error is exponentially small in as 
per eq. (|^. Naively, one might expect this to have a 
uniform distribution in position space, but Fig. shows 
that there remains a narrow ‘spike’ at the origin even in 
the high dimensional limit. This arises because loq only 
contains energy levels of odd n, all of which have maximal 
probability at the origin. 

Finally, armed with a time-evolving state Pcit) and 
its time-averaged state Wg we can investigate the equi¬ 
libration process. For the system and measurement in 
question, the time evolution of p), as defined in 

eq. (§, is displayed in Fig. |^for w = L/2. The effect 
of equilibration can be seen clearly for high dimension, 
as the system (i) starts sharply out of equilibrium, (ii) 
quickly drops to a small fraction of the initial distance, 
and (iii) oscillates within that distance for the predom¬ 
inant majority of the evolution. It is also remarkable 
that equilibration happens with such ease, where only 
extremely simple systems (deff < 50) are well distinguish¬ 
able from their equilibrium states. 

One common pattern in the literature [HISIIIIII] is to 
analyse how the infinite-time average of the distinguisha¬ 
bility scales with the effective dimension of the initial 
state. For our conditions. Fig. [^illustrates that it scales 
inverse-polynomially, at roughly l/2dgg, slightly faster 
than the known upper bounds for general systems [5], 
which scale with . This is likely because these 
bounds always involve a step employing Jensen’s inequal- 



FIG. 3. Log-scale plot of the time average of the distinguisha¬ 
bility for the Gaussian initial condition, (pG(t),‘^G), as 


per eqs. and ( |14[ ), for several values of deff- The line 

shows the best linear ht over the logarithmic data. 


ity, {'/Tfl) < ^ (d?^), which probably underestimates the 
scaling. 


B. The Time Scale 


From the evolution of DM„ipG{t),^G) it is easy to see 
exactly when the high-dimensional system equilibrates, 
given that after the initial decay it displays only a few 
short-lived fluctuations away from equilibrium. The red 
dots of Fig. [^ (the same data as the blue line of Fig. j^for 
a shorter time interval) display the equilibration profile 
quite clearly. Thus, for this system, the equilibration 
time scale can be defined as the characteristic width of 
the first peak. 

This leads us to our most important result of this sec¬ 
tion. As detailed in Appendix [^ in the limit deff ^ 1 
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FIG. 4. Numerical computation Dml (p(^)i‘^) compared to 

the leading order (eq. fH) and next-to-leading order approx¬ 
imations, for deft ~ 53. The expression for the next-to-leading 
order, as well as the full expansion, is found in eq. (A121 of Ap¬ 
pendix]^ The red dotted line is a zoom of the plot in Fig.]^ 
The constant region up to t/Tg « 3 x 10“^ happens because 
the wave function has not yet spread beyond the edges of the 
measurement window. 


and t Tcdcff, eq. ^ reduces to 


Dm I, (wg)Pg(^)) — 

2 TT 


where 


"6+0 



'G . - 


Tg = 


mL(j 
hiT ’ 


(15) 


(16) 


is the characteristic time scale for equilibration in this 
system. The significance of this expression is clear from 
Fig. a which compares the same numerical results pre¬ 
sented in Fig. a (for the high dimension case) with the 
expression above, in the short time limit where equilibra¬ 
tion occurs. For t « 0 you need higher order terms to get 
a precise estimate. However, for t « Tq, the lowest order 
approximation becomes the dominant contribution, and 
it is all one needs to describe the distinguishability in the 
vicinity of the equilibration time scale. 


One way of writing eq. (161 adimensionally is to com¬ 


pare it to the period of the groundstate Tg 


Tg 


1 


Tg Ih^TTfieff ’ 


(17) 


where Tg = = 


2'Kh AmL^ 




It is also interesting to compare this with recent re- 
suits [H] , which show that the equilibration time scales for 
typical measurements (random pojectors on the Hilbert 
space) scale with H/SDe, where SDe is the standard de¬ 
viation in energy of the state. In Appendix |B| we show 
that this corresponds to 


^typical _ I 


(18) 


FIG. 5. Numerical computation of the distinguishability, 
DuipJo, Poit)), averaged over time, for a range of different 
measurement widths. The line is the average value of this dis¬ 
tribution. Initial condition is a Gaussian wave-packet eq. Q 
with deff ~ 53. 


for our system. This proves that, while our system 
presents fast equilibration for projective position mea¬ 
surements, it is not nearly as fast as the time scales pre¬ 
dicted for random projective measurements. This corrob¬ 
orates the notion that measurements of physical signifi¬ 
cance (like the ones studied here) are much slower than 
the average. 


III. GENERALITY 


A. On the Measurements 

All the plots in the previous Section take w to be 
L/2. However, numerical tests were performed with var¬ 
ious different widths and neither the overall behaviour of 
p) nor the equilibration time scale varied signif¬ 
icantly. As seen in Fig. the amplitude of 73 m„ only 
shows significant change when w is very close to 0 or L — 
in which case it falls to 0 as the measurement conveys 
no information. This is evidence that these results are 
not a consequence of some particular symmetry emerg¬ 
ing from the chosen width, but are the behaviour one can 
expect from projective measurements in position space. 
In fact, most measurement of this class are slightly better 
than the w = Lj2 measurement, but this choice of w was 
helpful for our analytical results. 


B. On The Initial Condition 


In order to verify that our results are not a consequence 
of some special property of the initial condition, the same 
Hamiltonian and measurement were also investigated un¬ 
der the initial condition of a uniform distribution in en- 
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ergy space 


1^1^^) = 


VdcS 


deff 

E19 ■ 


(19) 


Satisfactorily, this system (described in more detail in 
Appendix presented a very similar equilibration pro¬ 
file. Its time scale was analytically estimated to be simi¬ 
lar to the Gaussian scenario in its proportionality to d~g, 


Tg 4(ieff ’ 


( 20 ) 


and the infinite-time average of its distinguishability 
scales with 0.77d^'^^, which is also very similar to the 
Gaussian’s scaling although the multiplicative constant 
is larger. 


IV. DISCUSSION 

Recent developments in the area have shown that 
generic time scale bounds are too weak to be physi¬ 
cally meaningful, while typical (but non-physical) mea¬ 
surements or typical Hamiltonians equilibrate extremely 
rapidly. As such, in the search for physical equilibration 
time scales, we must investigate the specifics of particu¬ 
lar systems and measurements. This work takes as step 
in that direction, thoroughly analysing a physical system 
and showing that fast time scales emerge, albeit much 
slower than those predicted for a typical observable. 


Here, we have exposed the equilibration dynamics of 
a wave-packet constrained inside a box when inspected 
by coarse-grained position measurements (onto a spatial 
window). We have derived an equation for the distin¬ 
guishability in the vicinity of the equilibration time scale, 
which also results in an equation for the time scale itself. 
This, in turn, is used to show that it scales with l/d^s, 
in constrast with the much faster time of 1/ d^g which is 
known for typical measurements in this system. 

Note that an estimate of the average time for the par¬ 
ticle to reach the edges of the box is given by 


(L/2) mLa 
" (Ap/m) " ~ir = 


( 21 ) 


Hence, perhaps surprisingly, the time scale for equilibra¬ 
tion is comparable to the time scale to reach the walls for 
the first time. In particular, this means that equilibra¬ 
tion is fast in a natural sense, and the particle does not 
need to ‘bounce back and forth’ many times in the box 
in order to equilibrate well. 

It is remarkable that, despite being a simple and well- 
studied textbook system, there are still questions to be 
asked about the particle in a box and still lessons to be 
learned. In particular, our results show that this sim¬ 
ple physical system can display equilibration behaviour 
over a physically realistic time scale, which reinforces re¬ 
cent developments concerning equilibration of abstract 
systems. 
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Appendix A: Equilibration Time Scale 

In this Appendix we use the following abbreviations, 
4> = 27rcr/L, Vnj = {En - Ej)/h, V = Ei/h, and pnj = 
{n\Po\j)- 

Claim 1. According to the set of outcomes M = 
{Ay^^ Byj\, with w = , the time-dependent distance 

between the state in eq. ([^ and its equilibrium state 
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(eq. 


2 (f> 

DM[u)G,PG{t)) « 

TT V27r 


kl 


cos(4fcZj^i)e-('="+'")TS 

kj — l 

(Al) 


sin(/c7r/2) sin(/7r/2) 

k I 


where B^i = [cos hr — cos kn] 

Where the approximate relation is exponentially precise 
in l/cj), as per eq. 

Proof. First, as defined in eq. © 


DMiuJc^Pait)) = Tr ALipcit) - ujg) 


{n\AL\j) 


(A2) 


= 2 


oo n—1 

EE Pn,j cos{vnjt) (n|Ai |j) 

n=2i=l 


where 
Pn,j ^ 


A(h pi P+rP niT jiT 

, _ e ^ « sin — sin — 

2 2 


(A3) 


7r(n — i) 7r(n + 7 ) 

cos-cos- 

2 2 


„ „ 31- 

2 


2 

ttJ. 

ii: 


sin {niTx / L) sin (jttx/L) da; 


sin (na) sin (ja) da 


cos ((n — j)a) — cos ((n + j)a) da 


sin {n — j) sin {n + j) 


n-j 


n + j 


(_l)n-a + i , f n-j 
-sin TT 


and note that 

{k + I + 1 )tt . {k — 1 )tt 


sin ■ 


sin ■ 


= 0 , 


yk,leZ, (A6) 


and 


{k + l\AL\k-l) = 


sin(7rfc/2) sin(7rZ/2) 

k I 


• (A7) 


Combining all of these identities with eq. (A2), yields 

D A;—1 


4 0 


TT 


Y Y {^klnt)e-^'^^+^"'>'^Bki 

fc=2 Z=1 


(AS) 

which can be easily reduced to eq. (|A1 ) by noting that 


the summand is symmetric with respect to k and I, and 
that Bkk = 0. □ 


Next, note that the double sum in eq. (All can be 


divided into two expressions, respectively proportional 
to sin(Z7r/2) and to sin(fc7r/2). Furthermore, sin(Z7r/2) 
is non-zero only when I is odd, which means (cos hr — 
cos/ctt) = —(1 -|- cos/ctt) which, in turn, means that the 
only non-zero terms of this expression those with odd I 
and even k. Then, the part of the double sum in eq. (A 81 
which is proportional to sin(Z7r/2) looks like 


-2EE cos {Aklvt)e 


-(k^+P)^ 


Furthermore, it is straightforward to see that 

(nlAjlj) = I If sin (n^(f + i)^ sin di 


/c=2 Z=1 
even odd 


sin(l7r/2) 

I 


■ (A9) 


Applying the same logic to the expression proportional 
to sin(fc7r/2) yields 


00 00 


-2EE cos {Aklnt)' 


-(k^+P)^ 


fc=i 1=2 

odd even 


sin(fc7r/2) 


, (AlO) 


which is identical to eq. ( |A9[ ) up to a swap of dummy in¬ 
dices. Thus, by writing the odd index in each expression 
as 2p -(- 1 and the even index as 2o, we have 


A’m(wg,/Og) 


(All) 




Y cos (8o(2p -k l)iyt) 


0=1 

P—0 


Sin TT 


(A4) 


2p+l ^ ’ 


7c(n-j) V 4 
(_l)n-Hj f_n + j 

TT{n + j) 

Finally, we state the identity 

OO n —1 00 fc —1 

EE= EE f{k + l,k — l)+ Dm(p^G)Pg) 

n—2j—1 k—21—1 

00 k — 1 

YYfik + l + hk-l), (AS) 

fc=lz=o 


,-(4o2 + (2p+l)2)^ 


We can further abbreviate this by defining i?„ = 


7r-\/^ 

2 Af 


cos (8(2p -I- l)ovt)Rpe 


0=1 

p—0 


V2tt 


OO 00 


'Y.RpYp cos (8(2p -I- \)ovt) 

p—0 0—1 
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which then produces eq. In particular, note that 

|i?i| = |i?o|®/3, which is why the p = 0 term is sufficient 
for a very good estimate of the equilibration time. On the 
other hand, at t = 0 this is very close to an alternating 
harmonic series, which is why it converges slowly and 
needs many terms to describe the constant region from 
Fig. El 


Appendix B: Energy Standard Deviation 

In this Appendix, we calculate the energy standard 
deviation (named here ve) for the Gaussian initial state. 
Taking ^ = 2^, we have, for any p e N, 


Tr[i7Vo] = EHPo|n)£;^ 


n—1 

2(1) 


sin^ ^ 


mT\ 

T ) 


= El^ y n2Pe-"'^sin2f —) 

n— — oo 

= Sp{V2(^), (Bl) 


where 


^ E (B2) 

y/TT 


]=-00 


Given that 


1 

Sp{l) = ^ dx x^Pe-"' + O y) , (B3) 

V ^ J — oo 

one has that, for small 7 , 


^1(7) = 2 + 0(7") 
^2(7) = ^ + 0(7"). 


(B4) 


FIG. 6. The wave function as defined in dClll, for different 
values of deff. 



X 

L 

FIG. 7. The spatial distribution of the equilibrium state Wu 


for the initial state in eq. (Cl I, for different values of deB- 


And finally, 

z;| = Tr[i77] -(Tr[i7p])' 

^ ^ oy 

4 - 


AEl 


4 22 

2 


= 4(+A') (TlI 1+0(0=) 

\2mL^ J \2 tt(7 J 2 ^ 


Ve = 


1 1 

2 w? 16(t"^ 


oy) 


4:^/2ma‘^ 


-id- 


(B5) 


Dividing H/ve by Tg = 7^ yields eq. (|l 8 l. 


Appendix C: Uniform Distribution 

The second initial condition considered was the con¬ 
stant distribution over N energy eigenstates, with 

1 ^ 

l'*^-) = ^El*)’ Pv = \d'.){y, (Cl) 


y/Ny 





































FIG. 8. pv{t)) with the uniform initial condition for 

two different defif- 
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FIG. 10. Log-scale plot of the time average of the distin- 
guishability for the uniform initial condition, DM(pu(t),<^u), 
as per eqs. ([^, (|G1[) and (G3l for several values of deff. 



FIG. 9. Short time comparison between the numerical com¬ 
putation of Dm{s^v, Pv{t)) (for deff = 500) and the g{t) given 
in eq. ([^. 


for some S N. This initial condition was chosen be¬ 
cause it allows for the exact evaluation of the most rele¬ 
vant part of the sum which defines the distinguishability. 
It is trivial to verify that deff = N. Fig.j^shows the prob¬ 
ability distribution of this state in position space, while 
Fig. [^displays the same for its corresponding equilibrium 
state ujv = jj En=i 

Since the initial state is clearly concentrated on the left 
side of the box, we use a slightly different measurement 
M = {ni,n/j}, which checks on which side of the box 
the particle is 

Ul= [ |a:)(x| dx, Un = 1 - n^. (C2) 

J -^/2 


2 fc 

with 1 / = 2 mL^ ■ This function is displayed on both curves 
of Fig. and on the blue curve of Fig. The equili¬ 
bration behaviour is, again, clearly visible in the high¬ 
dimensional limit, and we see (from Fig. that the first 
zero of the Dm function is a good approximation for the 
equilibration time scale. 


Following similar calculations as those from Ap¬ 
pendix]^ one finds that, for short times t <^Tg, eq. (C3l 
can be well approximated by the function 


9{t) 


1 

Ntt 


sin{2Nut) 

sin(^f) 


cos(:^f) 


(C4) 


Fig.|9] shows how these two functions compare to each 
other. To estimate the first zero of DM{t), we calculate 
the first zero of ^(t), which is exactly at 2sin(27Vi/t) = 
sm{2i/t). For large N, this is approximately 


sm{2Niyt) « i/t 
t ~ 

ANv AN' 


which yields the time scale 

Tu _ 1 ^ 1 

Tg AN 4deff ’ 


(C5) 


(C6) 


Again, we wish to study the time evolution of the dis¬ 
tinguishability according to this measurement, it is now 
given by the sum 


T>m(wu,Pu) = 




E 


n^j=l 


N 


[n\T^L\j) 


(C3) 


in good accordance with the time scale results for the 
Gaussian. 

For completeness. Fig. [T0| shows that the infinite-time 
average of the distinguishability in this scenario scales 
with d~g' , which is very similar to the Gaussian’s 
although the multiplicative constant is larger. 








































